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I. 



INTRODUCTION 



Collapse models are typically described by stochastic differential equations (SDEs) which involve 
white noises at infinite temperature. Three important examples are given by the GRW model [H[2]> 
the first collapse model, the CSL model [21 [3], which is the generalization to field theory of the 
previous model, and the QMUPL model [HE]. This last model has the advantage of being physically 
reasonable and at the same time mathematically analyzable in detail [6l414j. 

If one wants to identify the noise causing the collapse with a physical field, it is necessary to 
generalize the previous models. The reason is twofold: first, a white noise has a flat spectrum, i.e. 
every frequency appears with the same weight, while for a physical field one would expect to have 
a colored spectrum, possibly with a high-frequency cutoff. Secondly, the infinite temperature of 
the noise is clearly an unphysical feature. 

The first of these two issues is solved by non-Mar kovian collapse models |15H18j . These models 
involve Gaussian noises which have a general time correlation function, i.e. a colored spectrum. 
Concerning the second issue, a generalization of the QMUPL model to finite temperature noises 
have been studied in [19]. At the practical level, this goal is achieved introducing in the SDE new 
terms which account for dissipation in the interaction between the system and the (still white) 
noise. So far, these two generalizations have been studied independently, because of the difficulties 
in the calculations involved. 

Aim of this paper is to study the generalization of the QMUPL model both to a non-white and 
finite temperature Gaussian noise. Such noise displays, at the same time, both the features which 
make it comparable to a physical field. This is then the ultimate generalization for this model, 
and it represents the first description of a collapse model by means of a noise which has physical 
features. 

The SDE describing the evolution of a wave function according to the non-Markovian dissipative 
QMUPL model is readily obtained substituting in Eq. (11) of |20j the operator L = q + i^p, where 
q and p are the position and momentum operators respectively. The positive constant /x accounts 
for dissipation, and it is linked to the temperature T of the noise by the relation [19] : 



ti 2 



(1) 



T = 



4:ink B fj, 



The equation for one particle reads: 
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where Hq is the Hamiltonian, {•,•} denotes the anticommutator, A is the collapse parameter and 
w(t) is a non-white Gaussian noise whose time correlation function is D(t,s). The form of the 
operator L, and the introduction of the anticommutator, are justified by the fact that we want to 
reproduce known results in the white noise limit |19| . The generalization of this equation to many 
particles systems can be easily derived |17j . One can also check that in the limit jj, — > one obtains 
the SDE describing the non-Markovian QMUPL model \W\ I17j. and in the limit D(t, s) — > 6(t — s) 
one recovers the dissipative QMUPL model [19]. When both limits are taken, one recovers the 
original QMUPL model [3]. 

Goal of this paper is to give the analytic solution of the SDE ^ for the harmonic oscillator, 
and to analyze the collapse process. We remark that, like in all collapse models, Eq. ^ does not 
preserve the norm [2]. In order to obtain the physical states ipt ° ne nas to normalize the solution 
<Pt of Eq. (§: 



^=niV ( 3 ) 



The physical probability distribution is determined by a measure P, which is defined as follows [21 

dP = ||^|| 2 dQ, (4) 

where Q is the measure according to which the properties of the noise w{t) are defined. 

This paper is organized as follows: in section [IT] we compute the Action associated to Eq. Q; 
this Action will be used in the path integration through which we will find the solution of Eq. ^ . In 
section [TIT] we explicitly evaluate the path integration and we write the Green's function associated 
to Eq. ([2]). In order to better understand the properties of the dynamics, in section IV we consider 
an exponential correlation function D(t, s), and we analyze the behavior of Gaussian wave functions. 



II. NON-LOCAL ACTION OF THE MODEL 

In order to solve Eq. ^ we will proceed as first suggested in [201 EI] arid done in [161 E! : we 
will find the associated Green's function G(x,t; xo,0), 

/oo 
dx G(x,t;x ,0)4> (x ) , (5) 
-00 

and we will derive the explicit expression for the propagator. In [21] it has been shown that it is 

possible to write the Green's function via the path-integral formalism: 

rq(t)=x rq(t)=x t 

G(x,t;x ,0) = / V[q]e s[q] = / V[q]V\p] el Jo ds P (s) q (s)-H(s) ^ (g) 

Jq(0)=x Jq(0)=x 
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where T>[q] represents the functional integral over all the paths connecting q(0) = xq to q(t) = x, 
and the functional S[q] is the classical Action associated to the system. The right hand side of 
Eq. Q refers to the Hamiltonian formulation of the path- integral formalism |22j . where q is the 
velocity and H(t) is the classical Hamiltonian associated to the system. One can easily check that 
substituting in Eq. ^ the following Hamiltonian 

H(t) = H + ^{q,p} + ihV\(q+^p\w(t)-2ih\q^ dsD(t, s) (q(s) + , (7) 

one obtains the correct propagator associated to Eq. ([2]). Since from a technical point of view it 
is easier to compute the Lagrangian path integration of Eq. ^ instead of the Hamiltonian path 
integration, we have to determine the Action 

S[q] = f dsp(s)q(s)-H(s), (8) 
Jo 

starting from H(t) of Eq. ([7]) and transforming the dependence on (q,p) in a dependence on 
(q,q). One can try to apply the standard Hamilton formalism to Eq. ([7]), in order to compute 
the Lagrangian from H(t), but soon realizes that a problem lies in the fact that H{t) is a time 
non-local Hamiltonian, i.e. it contains a memory term which accounts for the whole past history 
of the system. In particular, in this case it becomes problematic to evaluate 

for the integral term of Eq. ([7]). Since in local Hamiltonians only the time t appears, there is not 
ambiguity in taking this derivative. One can then argue that only the end-point t of the integral 
in Eq. ([7]) contributes in Eq. ([9]). However, also the Action computed with this approach does 
not give the correct result. In the essence, when dealing with time-non-local Hamiltonians, the 
standard formalism does not apply anymore, and one has to resort to a new formalism. We then 
had to setup a new formalism based on variational calculus |23j . Such formalism is a generalization 
of the standard Legendre formalism, in which one makes use of the variational calculus to write 
the equations of motion. Basically, the derivatives with respect to the position and momentum 
coordinates, are substituted with functional derivatives of position and momentum at a given time. 
In this formalism, for example, the generalized Euler-Lagrange equations read: 

5S[q, q] d 5S[q, g] 
5q(s) ds 5q{s) ' 1 ' 

note that the presence of the Action S[q, q] instead of the Lagrangian is crucial because of the 
possibility to get some infinities from the functional differentiation. Moreover, like in the Euler- 
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Lagrange formalism, q and q are independent variables. Furthermore, introducing 

r-t 



n[q, P ] 



H(s) ds 



o 



one generalizes the definition of conjugated momentum as follows: 

sn[q, P ] 



?00 



8p(s) 



(11) 



(12) 



In this way, one can conveniently write the non-local Action S[q, q] associated to the Hamilto- 
man 0. Note that this is a non-standard Action, since it involves imaginary terms other than a 
double integral term. These terms are eventually responsible for the collapse of the wave function. 
Substituting Eq. Q and Eq. ^ in Eq. after some calculations one finds that for an harmonic 
oscillator, S[q,q] takes the following expression: 

ds- — q 2 (s) — mXuq(s) q(s) — 



S[q, 



H 



rn 
~2 



q 2 (s) — mXjj, q(s) q(s) — — ft 2 q 2 (s) + mVX/i q(s) 



ihVXw(s) + mX 3/2 fi 2 w(s) +2mA 3/ V 2 / drD(r, s)w(r)) q(s 



(2mA 2 /i 2 + 2ihX) q(s) / dr D(s,r)q(r) — 2mX^q(s) / drD(s,r)q(r) 
Jo Jo 



+2mX /i / dr D(s,r)q(r) / dr' D(s,r')q(r') + 



mXfi z 



w 2 (s) 



(13) 



with U 2 = co 2 — X 2 [i 2 , where u is the frequency of the oscillator. This is the action corresponding 
to the Hamiltonian of Eq. ([7]) . 

III. SOLUTION OF THE EQUATION FOR A QUANTUM HARMONIC OSCILLATOR 



Aim of this section is to compute the path integration of Eq. ^ with the Action given by 



Eq. (13). The main result of this section will be the Green's function given by Eqs. (44)-(49). 



First of all we rearrange the Action of Eq. (13), collecting terms of the same kind: 

ft 



where 



ds- 
o " 



m 



2 Q 2 ( s ) ~ mXfJ,q(s) q(s) — —Q 2 q 2 (s) + mVXfiq(s) 



/•S r-S 

— A(s) q(s) + q(s) / dr B(r,s)q(r) — 2mXfiq(s) / drD(s,r)q(r 
Jo Jo 

A(s) = ihVXw(s) + mX i/2 n 2 w(s) + 2mX 3/2 fi 2 [ drD(r, s)w(r) , 

Jo 

rr 

B(r, s) = (2mA V + 2mA) D(r, s) + 4mA V / dr'D(r, r')D(s, r') . 

Jo 



(14) 

(15) 
(16) 
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We remind that in the path-integral formalism q and q are not independent and therefore the Action 
is a functional of q only. We also stress that, in spite of many similarities, there are important 
differences with respect to the calculation we performed in the dissipation- free case |16| [T7] : terms 
depending on q(s) appear both coupled to the noise w(s) and inside the double integral. Secondly, 
the double integral terms range from zero to s (instead of t), and the function B(r,s) is not 
symmetric. All these features of the Action make the path integral much more difficult to solve 
than that of [T2l[T6]. 

We will proceed following the standard polygonal approach proposed by Feynman \22\ l24l [25] . 
The time interval [0,i] is divided in iV subintervals, each of length e = t/N, and the intermediate 
time points are naturally defined as t k = ke. One can then define a "discretized propagator" as 
follows: 

N , ,.N-1 

G N (x,t;x ,0)=(^-) 2 ■■■ ]Jdq k e s ^, (17) 
m J J k=i 

where S]\r[q] is the discretized form of the Action S[q], and q k = q(t k ), k = 1, . . . , N — 1. The path 
integral is then given by the limit N — > oo of the multiple integral over the N — 1 variables q k and, 
as a consequence, 

G(x,t; xo,0) = lim Gjy(x, t; xo, 0) . (18) 

N— >OD 



Unlike the non-dissipative case, the Action of Eq. (14) has linear terms depending on q. It has 
been shown |22j that in this case the discretized Action cannot be evaluated, as it is done in 
the standard case, at the initial point of each time subinterval, but it must be evaluated at the 
mid-point Xh . This is the so called "mid-point formulation" , and in this case the discretized 
Action reads: 

N 

h 



S N [q] = £ 1 
k=l 



m. .o m\u, w , rn o 2 

— {q k - q k -i) — {qk + qk-i){qk ~ qk-i) ~ e— S2 (q k + <j fe -i) 

Z€ Z o 



A + 
+mV\/jw k (q k - q k -i) - e^(q k + q k -i) + e 2 <lk j^ ^Bj^qj + q^i) 



k 

-em\n(q k + q k -i) ^ D jtk (qj + g,_i) 



(19) 



where w k = w(t k ), A k = A(t k ) and Bj ;k = B(tj,t k ), Dj )k = D(tj,t k ). The constraints of the path 
integration are: qo = xq and q^ = x. The most convenient way to perform the iV — 1 integrals of 



Eq. (17) is to reduce them to Gaussian integrals, which can be easily solved. It is then useful to 
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write Gn(x, t; Xq, 0) as follows: 

N 



m mX/j, m^ 2 ^- 

— e—n +C NjN x +Cl N x : 



2c 



m mXfj, m 



+ [ 2c + 2 



e-0 2 + C+ )x 2 + K^x + K+: 



/+oo 
dX. exp [-X • MX + 2 X • Y] , (20) 
-oo 

where we have introduced the following vector notation (T denotes the transpose): 



X = (q 1 . . .q N -i) , 



(21) 



' / iii in 2 Ci,i 



with 



h \ 2e 



+ e—n z 



+ 



i 1 
H2 



- \ 



K+ + K^ 



V K N-1+ K N ) 



(x Q \ 



\ / 



+ 



ft 2 



ft I 2e 



" v~» 7 



( c+ A + c+ ^ 

£%2 + ^3 

^l,N-2 + ^1,N-1 
V ^1,N-1 + ^N-1,N / 







I 
n 



C 



N-2,N 



A? 



,(22) 



V ^N-1,N + ^JV.AT / 



J 2 



2 



(23) 
(24) 



The matrix M is the sum of two (N — l)-dimensional square matrices M and M, whose entries 



are: 



Mi 



Mi 



i,i±l 



c j+l,j) 



2e 8 



2 3,0 



M id = , j / i, i ± 1 , 



(25) 

(26) 
(27) 



i<j = ~^2 + ^i+M+i + ^J'.i+i + ^iW) ' 



(28) 
(29) 



s 



These two matrices have this nontrivial expression because of the features of A(s) and B(r, s) of 



Eqs. (15)-(16), previously described. The integration of Eq. (17) is evaluated using the standard 



formula for multiple Gaussian integrals: 

30 

dX exp [-X • AfX + aX • Y 



7T 



N-l 



y det(M) 

The discretized propagator Gn(x, t; xq, 0) then becomes: 



cxp 



-Y • M~ l Y 



(30) 



G N (x,t; x ,0) 



N 



'(m/2itve) 
vrdet(M) 



exp 



Y • M~ l Y + 



i ( m mXfi 



h V 2e 



i I ' m mXn 



h\2e 



+~K+x + ~KQ X0 + O(e) 



(31) 



where 0(e) collects all the terms of order e or higher. In order to have the exact expression for 
G(x, t; xq, 0) we need to compute Y • M~ 1 Y and take the limit for N — > oo. We introduce the 
vector Z := (zi, . . . Zjv_i) T , which is defined as follows: 



Y Z = Y M~ X Y 



(32) 



i rn m 2 _C^ 
h\2e 8 2 



i I m m 



XqZi 



N-l K + , K - 



h\2e 

N-l n + 



c 



N,N 



XZ N -l 



. N-l r - i ^»+ 



Zi 



j=l j=l j=l 

We can then assume that the components of Z are the discretized values of a function z(s): 



Zk = z(tk), k = 1. . .N — 1, and z(0) = xo, z(t) = x. Substituting this expression in Eq. (32) one 
finds: 



Y • M~ l Y 



-|^xo (*(0) + 6i(0) + 0(e 2 )) - l ^-x (z(t) - ez(t) + 0(e 2 )) 



2V-1 



AT-1 



3=1 3=1 



i emXu \— ^ 



AT-1 



ft 2 



(33) 



3=1 



where the dot denotes the derivative of z(s). Substituting Eq. (33) in Eq. (31) we obtain: 
G N (x,t;x ,0) = \l '_ exp 



vrdet(M) 



ira 



2h (x z(0)-xz(t)) +h 2 



1 ,,_ i '"-V AU . . 

+ -K x + t_ — - — 2^ ( W J - ™3+l) z 3 



i m 



\/XfJ> 



JV-l 



ft 2 



3=1 



ft 4 



N-l 



J2 ( A i + z i 



3=1 



+ 



i emX/i 
ft - 2~ 



N-l 



Xq 



J2(D ltj + D 1 , j+l )z j + 0(e) 



3=1 



(34) 
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The last step is to compute the determinant of M. Since the calculation of det(M) follows the 
same procedure as in |17l I26j . we will report only the result: 



(vn \ ^ — i 
me) UN ® 



(35) 



where UN(t) is a suitably defined function, which in the limit N — > oo converges to u(t). Note that 
u(t) enters the global factor: the real part is re-defined after one renormalizes the wave function in 
order to obtain physical states. The imaginary part gives an unimportant phase factor. Hence, we 
are not interested to the explicit expression for u{t). Substituting this result and taking the limit 
N — > oo, one finds that the Green's function takes the expression: 



G(x,t;x ,0) 



in 



cxp 



m 



xoi(O) — xz(t)) 



mXfi , 2 2 s 
[x - x ) 



+mv Xji \w{u)xq — w[t)x) H — 



ds w(s)z(s) 



o 



ft A(s) /"* 
+ / ds — - — z(s) — mA/iXo / dsD(0, s)z(s) 



o 







(36) 



where the function z(s) will be evaluated now. Considering the definition of Eq. (32), once written 
in components and divided by e, one obtains: 



in 



TO. 



mX/i 



2e2 {zk+i ~ 2z k + z k -i) + — tt (z k+1 + 2z k + z k -i) + -^-D kk (z fc _i + z k+1 ) 



mX/j, 



mX/j, 



-D kk + D k+ \ k+ i - D kk+ i - D k+lk ) z k - - ^ (B k j + B k+ ij + i + B k j +1 + B k+ ij) Zj 

j=i 

k N-l 

^ (Ac+ij+i — D k j + i + D k+ ij — D k j) zj — - ^2 (Bjk + Bj +lk+ i + Bj k+ i + Bj + i k ) Zj 



mX/j, 



3=1 
N-l 



j=k 



(Dk+ij+i — D k j + \ + D k+ \j — D k j) Zj 



j=k 



mVXfi w k+1 - w k A k+1 + A k mXfi 
= ^ ~ e 4 + ~~ 2~~ X ° ^ lk+1 > ' 

Taking the limit e — > one finds that z{s) satisfies the following equation: 



U - z(s) + | {n 2 + ±X»D(s, s)) z(s) -\fdr B(r, s)z{r) - \ f 

~ J J s 

f s dD(r,s) f 1 dD(r,s) 
-mX[i I dr — - — z(r) — mXfi / dr — - — z{r) = 



k = l,...N-l. (37) 



dr B(s, r)z(r) 



o 



mV~Xfi A(s) 

o — w { s ) n ^ mA A* x D(0, s) . 



(38) 



Because of the boundary conditions we imposed on z(s), the solution of Eq. (38) depends on 



the end points xq and x. In order to make this dependence explicit in the expression of the 
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propagator (36), we rewrite z(s) as follows: 

z(s) = f(s)x + g(s)x + h(s). (39) 
Defining the following integro-differential operator: 

I[e(s)} := — e(s) + — (O 2 + 4A/x.D(s, s)) e(s) - 1 [ dr B(r, s)e(r) - 1 [ dr B(s, r)e(r) 



-mA/i / <ir 



dD(r,s) 



dr 



e(r) — raA/i / dr 



&D(r, s) 



<9s 



e(r) , 



(40) 



one can prove that f(s) satisfies the following integro-differential equation: 

I[f( S )] = m\fiD(0,s), (41) 
with boundary conditions / (0) = 1, f(t) = 0, while g(s) solves the homogeneous equation 

l[g(s)) = 0, (42) 
with boundary conditions g(0) = 0, g(t) = 1. The function h(s) instead satisfies the non- 



homogenous equation: 



I[h(s) 



m 



■wis) 



A(s) 



(43) 



2 w 2 ' 

with boundary conditions h(0) = h(t) = 0. Note that /(s), ^(s) and h(s) depend also on the 
parameter t, so we will add the subscript t when confusion could arise. Here and in the following, 
the dot denotes differentiation with respect to the variable in parentheses. According to the ansatz 



we introduced, we can rearrange the propagator (36) as follows, making the dependence on xq and 
x explicit: 

A t xl - A t x 2 + B t x x + C t x + V t x + £ t 



G(x, t; xq, 0) 



in 



exp 



(44) 



Imfit u(t) 

The coefficients At, At and Bt are deterministic functions of time, while the coefficients Ct, T>t and 
£t are random processes: 

At = k (f t (0) - A/i - 2A/i / ds D(0, s)f t (s)] , A t = -k {g t (t) - A//) 



Bt = k( f t (t) - g t (0) + 2A/i / ds D(0, s)g t (s) 



k 



im 
2K' 



C t = -k (ht(0) + 2V\fiw(0) + V^M^ 



+ 



A(s] 



V t = k (h t (t) + 2V\fiw(t) - \/Xfj, J 



m 
t 



w(s)ft(s)ds 

ft(s)ds - 2\fi [ dsD(0,s)h t (s] 
Jo 



w(s)g t (s)ds 



St. 



k vA/j / w(s)ht(s)ds + 



A(s) 



m 



Ms] 

rn 

t 



gt(s)ds 



h t (s)ds — Xjjt I w (s) ds . 



(45) 
(46) 

(47) 
(48) 
(49) 
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This is the result we wanted to arrive at: the explicit expression of the propagator associated 
to Eq. ([2]), for an harmonic oscillator. This is the main result of the paper. It is worthwhile 
noticing that the propagator has a Gaussian form like in the white noise case: this implies that 



Gaussian states evolve preserving their structure. We stress that under suitable limits Eqs. (41)- 



(49) reproduce known results. 



IV. EXPONENTIAL CORRELATION FUNCTION 



In the previous section we computed the Green's function of Eq. ([2]) for the quantum harmonic 



oscillator in terms of the solutions of the integro-differential equation (38). However, such an 



equation cannot be explicitly solved for any D(t,s). In this section we consider the case of a 
correlation function which allows to solve Eq. ( |38[ ) explicitly, and moreover it is also physically 
meaningful. This is the exponential function: 



D(t, s) =- - < 



l e -y\t-s\ 
2 



(50) 



where 7 is the frequency cutoff. First of all we note that 

dD(r,s) dD(r,s) 



ds 



dr 



(51) 



which allows to rewrite two of integral terms of Eq. (40) as follows: 

dD(r,s) 



-2A// r dr dD ^ S K {r) - 2A/x f dr 







dr 



ds 



-e(r) 



2Xfi (d(0, s) e(0) - 2D(s, s) e(s) + D(t, s) e(t) + J D{r, s)e(r)dr - J D(r, s)e(r)dr^J .(52) 



Substituting this expression and Eq. (50) in Eq. (40), using Eq. (41) one finds that f(s) satisfies 
the following integro-differential equation: 

'3A 2 /i 2 7 ih\^ rt 



f(s) + n z f( s 

ft 



+ 



-A^7 / ( 

' s 



m j jo 

- f(r)dr +^-L [ dre-^ r+s) f{r) =0, 



dre-^ s - r \f{r) + \^ / e -^ s ~ r) f{r)dr 
Jo 



(53) 



2 Jo 

The procedure used to find the solution of this equation is the following [28J : we differentiate 



Eq. (53) twice, obtaining: 



/ (s) + {Q 2 + 2A/i 7 ) f(s) + 3AV 2 7 2 + 



m 



f{s)+ 1 2 



-y(s-r) 



A/X7 / ( 

AV 2 7 " 



f(r)dr 



(54) 



'3AV7 + ^r}X 1 \ /" dr e - 7 | s _ H /(r) _ A f* e - 7(r - s) j (r)dr + 
. 2 m J Jo Js 2 



dre-^ r+s ^f(r) 
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and substituting Eq. (53) in it, one arrives at the following fourth order differential equation: 



'f'(s) + {n 2 + 2A/x 7 - 7 2 ) /(a) + ( 4AV 2 7 2 - 7^ + 
The general solution of this equation is 



m 



M = o. 



(55) 



f(s) = fi sinh v\s + f2 sinh V2S + cosh v \s + cosh V2S . 



(56) 



The roots of the characteristic polynomial associated to Eq. ( 55 ) are not easy to handle. For the 



sake of simplicity, we limit our analysis to the case of a free particle (u = 0). In such a case, the 
roots become: 



Vl,2 



[(7"A/i) 2 ±C], 



( = J ( 7 - A//) 4 - 16A 2 /x 2 7 2 



8ih\~f 2 



rn 



(57) 



Since this is a fourth order differential equation, we need four boundary conditions to determine the 
coefficients fi univocally: two of them are /(0) = 1 f(t) = 0, while the other two are determined 



as follows. We evaluate Eq. (53) in and t: 



f(0) - ( AV 2 7 + ^) f dl e^ l f(l) - A M 7 f dl e^ l f(l) = , 



/(*) 



+ f dle -^ f{l) (58) 

2 rn J J 

+A/i 7 f + fdle-^f(l) = 0, 



and we substitute in the integral terms the expression for /(/) obtained rearranging Eq. (55). 
Integrating by parts, after some calculations one finds that the boundary conditions are: 



/(0) = 1, 



/(*) = o ; 



, ■ > ■ -> , ih\\ 

A V + A/X7 + 

m J 



iHX 



f (0) + A M7 /(0) = 7 ( 2AV + A/U 7 + — ) /(0) + A"VV, (59) 



f(t) + 2X^f(t) = - 7 f(t). 



The equation for g(s) can be found in the same way using Eq. (42), and it results to be: 



3A> Z 7 i^A7 



-A/X7 



/ Jo 

v2,,2„, 



(60) 
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Using the procedure previously outlined, one can easily show that also g(s) satisfies Eq. (55). 



However, the coefficients <?, entering the solution are not the same as fi since g(s) satisfies different 
boundary conditions, namely: 

9(0) =0, <?(*) = 1, 



AV + A/i7 + —) [9(0) + 2\^g(0)} = 7 f 2*V + A^ 7 + — ) g(0) 
m J \ m J 



g(t) + 2\^g(t) = - 7 g(t) - 2A^ 7 2 



(61) 



Applying once again the same procedure to Eq. ( 43 ) , one can derive the fourth order differential 



equation satisfied by h(s), whose solution can be found generalizing the technique used in |17j . Since 
the analytic expression of h(s) is not necessary for the analysis of the next section, we omit the 
calculation. 



The solution of the systems of equations (59)- (61 ) are very long and it does not make any sense 



to report them here. However, their expressions can be easily derived by using standard software 
like Mathematica® . We stress once again that the final solution reproduces known results under 
suitable limits. 

We eventually stress that the expressions for f(s), g(s) and h(s) provide the explicit expression 
for the solution of Eq. Q. For the case of an exponential correlation function every coefficient is 
determined. This is a highly nontrivial result for a (infinite dimensional) non-Markovian dissipative 
collapse model. 

V. EVOLUTION OF A GAUSSIAN WAVE FUNCTION AND COLLAPSE FEATURES 



In order to analyze the collapse features of the model, in this section we study the evolution of 



Gaussian wave functions. Being the propagator (44) Gaussian as well, the shape of Gaussian wave 



packets in preserved during the evolution. The wave function at time t is given by: 

(p t (x) = exp[-a t x 2 + fax + 74] , 



where: 



a t = At 



4(a + A t 



Ct + Po 

2(0*0 + A) 



B t + V t , jt = 7o + £t + 



(Ct + Pof 
4(a + A t 



(62) 



(63) 



and the functions At - £t have been defined in Eqs. (45)-(49). One can see that the evolution of 



the spread of the wave function is deterministic, since At, At and Bt do not depend on the noise 
w(t). 
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FIG. 1. Time evolution of a{t) for a 1 kg particle with initial spread Co = 10~ 7 m, for different values of /i 
(m 2 ). Smaller values of \x imply stronger collapse. On this scale, the plot for [i — coincides with that for 
/j = 1CP 32 . Time is measured in seconds, distances in meters. 



We now focus on the case of an exponential correlation function, since in this case we can 



explicitly determine all the coefficients of Eq. (62 ). In particular, we are interested in the evolution 



of the spread a{t) of the wave function, since this is the quantity which shows how the collapse 
mechanism works: 

a(t) = , (64) 

where the apex r denotes the real part. The analytic expression of a(t) can be derived from 



Eqs. ri63h, rt45| and (56). 



Figure [T] shows the time evolution of the spread of the wave function of a 1 kg particle with 
initial spread gq = 10~ 7 m, for different values of /i. First of all, one can see that the wave function 
shrinks in time, and the behavior of the spread is qualitatively the same as in the non-Markovian 
model (infinite temperature) [TBI H7] . Moreover, as expected, the stronger the dissipation, the 
lower the temperature of the noise and the weaker the collapse process. Fig. [2] shows the behavior 
of a(t) at a fixed time t = 2 x 10 5 s, for different values of 7. As we could expect, as a general 
behavior the larger the value of 7, the more frequency components of the noise enter into play 
and the faster the collapse of the wave function. However, one can see that the a{t) is not strictly 
decreasing with 7. Moreover, the comparison of the dissipative non-Markovian model (solid line) 
with the non-Markovian one (dashed line) shows that dissipation makes this transition smoother. 
This is due to the fact that dissipation makes the interaction with the noise less effective (as one 
can also see in Fig. [I]). 
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FIG. 2. Dependence of a(t) on 7, for a fixed time t — 2 x 10 5 s. Initial values are m = 1 kg, er = 1CT 7 m, 
yu = 5 x 10~ 19 m 2 . The bigger 7 the stronger the collapse, unless for some values of 7. Time is measured in 
seconds, distances in meters. 



VI. CONCLUSIONS 



We investigated the dynamics of an harmonic oscillator according to the non-Markovian dis- 
sipative QMUPL model. This first time a model involving a noise with physical features (finite 
temperature and colored spectrum) is studied. We provided the explicit formula of the evolution 
in terms of the Green's function of the process. This is a highly nontrivial result. We studied 
the evolution of Gaussian wave functions in the free particle case, focusing on the behavior of 
their spread. We showed that the wave function collapses qualitatively with the same behavior 
as previous models. Moreover, we showed that both non-Markovian and dissipative effects make 
the collapse process less effective. The analysis presented can be considered an arrival point for 
the study of collapse models, since the model analyzed accounts both for memory and dissipative 
effects. 
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